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基于 飞 蛾 火焰 优化 算法 的 PMSM 分 数 PF* 阶 研究 ， 


张 学 典 ， 孙 俊 峰 ,， 秦 晓 飞 
(上 海 理工 大 学 光电 信息 与 计算 机 工程 学 院 ， 上 海 200093) 


摘 要 : 针对 PMSM( 永 磁 同 步 电机 ) 分 数 阶 P 信 控制 器 的 参数 整定 , 提出 了 一 种 频 域 法 和 飞 蛾 火焰 优化 算法 相 结 合 的 方 
法 。 以 永 磁 同 步 电 机 的 分 数 阶 模型 作为 被 挖 对象 。 首 先 使 用 频 域 法 ， 根 据 系统 的 相对 稳定 性 和 增益 变化 的 鲁 棒 性 等 条 
件 得 到 控制 器 的 方程 ; 然后 结合 该 方程 的 图 像 求解 出 分 数 阶 PI 控制 器 的 参数 ， 以 该 参数 为 中 心 位 置 ， 指 定 寻 优 的 范 
围 ， 进 而 用 火焰 优化 算法 对 其 周围 进行 寻 优 。 仿 真 和 实验 结果 表明 ， 该 分 数 阶 PI 控制 器 能 够 降低 系统 的 超 调 量 ， 提 高 
系统 的 快速 跟踪 能 力 和 抗 干扰 能 力 ， 具 有 较 强 的 便 棒 性 。 
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Research of fractional order PI* for PMSM based on moths fire optimization algorithm 


Zhang Xuedian', Sun Junfeng, Qin Xiaofei 
(School of Optical-Electrical & Computer Engineering, University of Shanghai for Science & Technology, Shanghai 200093, 
China) 


Abstract: Aiming at parameter tuning of fractional order PI* controller of PMSM, the author proposes a method, which 


combines the frequency-domain method and moths fire optimization algorithm. This paper used the fractional order model of 


PMSM as the controlled object. Firstly, it used the frequency-domain method, the it obtained the equation of the controller 
according to the relative stability of the system and the robustness of the gain variation. After that, it solved the parameters of 
fractional order PI controller by combining the image of the equation and optimized these parameters by using moths fire 
optimization algorithm. The simulation and experiment results show that the fractional order controller can reduce the overshoot 
of the system, improve the rapidity and anti-interference ability of the system, and has a certain degree of robustness. 
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等 人 采用 Outaloup 递 推 滤波 器 近似 ， 开 发 了 CRONE 控制 器 ， 
且 证 明了 CRONE 控制 器 比 传 统 的 整数 阶 PID 控制 器 具有 很 

水 磁 同 步 电 机 已 经 广泛 得 应 用 到 了 人 们 生活 的 各 个 领域 。 大 的 优势 四。Matignon 等 人 针对 分 数 阶 系统 的 稳定 性 和 可 控 琐 
随 着 我 国 “ 中 国 制 造 2025” 制 造 业 强 国 战略 的 提出 ， 像 数控 机 “性 等 方面 进行 了 研究 站。1994 年 ，Podlubny 在 研究 分 数 阶 控 币 
床 、 机 器 人 等 领域 都 需要 电机 能 够 被 精确 的 控制 ， 而 永 磁 同 步 ”系统 基础 上 提出 了 分 数 阶 PPDr* 控制 器 外。 
电机 因 其 优异 的 性 能 被 应 用 到 了 这 些 领 域 作 为 运动 控制 的 关键 近年 来 分 数 阶 PPDr 控制 器 得 到 了 许多 学 者 的 关注 ， 胡 海 
部 件 053。 但 是 受 电机 参数 变化 、 负 载 扰动 等 不 确定 因素 的 影响 ， ” 波 等 人 外 结合 细菌 砚 食 算法 、 粒 子 群 算法 对 分 数 阶 PID 控制 器 
控制 系统 的 动态 响应 和 抗 干扰 能 力 仍然 不 能 满足 现代 工业 控制 ”进行 了 设计 ， 提 高 了 传统 粒子 群 算法 的 计算 效率 ， 同 时 得 到 的 
的 要 求 。 因 此 ， 如 何 保证 系统 具有 良好 的 动静 态 响应 性 能 和 较 分数 阶 PID 控制 器 比 整 数 阶 控制 器 具有 更 好 的 动态 特性 。 吴 振 
强 的 鲁 棒 性 ， 一 直 是 永 磁 同步 电机 控制 界 研究 的 热点 中 。 随 着 。” 宇 等 人 09 把 分 数 阶 PID 应 用 到 了 智能 车 控制 ,最 后 证 明 分 数 阶 
现代 控制 理论 的 不 断 发 展 ， 许 多 先进 的 控制 方法 被 提出 ， 特 别 ”PID 控制 器 具有 比 传统 PID 控制 器 更 好 的 动态 性 能 和 更 强 的 鲁 
是 分 数 阶 理论 的 出 现 ， 分 数 阶 控制 器 是 借用 整数 阶 逼 近 的 方式 。” 棒 性 。 秦 昌 诚 等 人 [利用 自 抗 扰 分 数 阶 PID 控制 方法 , 研究 了 
在 计算 机 中 实现 ,因此 很 多 学 者 对 其 进行 了 研究 ，Oustaloup ”高 超声 速 飞行 条 件 下 的 飞行 器 再 入 滑翔 段 制 导轨 迹 跟 踪 ， 表 明 
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了 分 数 阶 PID 具有 更 强 的 鲁 棒 性 。Altintas 等 人 03 结 合 遗传 算 
法 和 分 数 阶 PID 优化 了 磁悬浮 系统 控制 器 的 参数 ， 通 过 分 析 验 
证 了 FOPID 具有 更 好 的 灵活 性 , 使 系统 具有 很 好 的 动态 响应 特 
性 021。 

本 文 是 在 频 域 法 的 基础 上 ， 应 用 飞 蛾 火焰 03.49 优 化 算法 对 
分 数 阶 PT 控制 器 三 个 参数 进行 优化 。P 控制 器 的 典型 公式 如 
式 (1) 所 示 。 


C(s) = 有 0+ 多) (1) 


其 中 : ,是 分 数 阶 比例 和 积分 参数 ，4 为 分 数 阶 阶 次 。 

飞 峨 火焰 优化 算法 四 是 2015 年 Mirjalili 等 人 提出 的 一 种 新 
型 仿生 群体 智能 算法 ， 与 其 他 的 智能 优化 算法 相 比 ， 它 具有 搜 
索 速 度 较 快 、 调 节 参 数 少 、 收 敛 精度 高 和 和 鲁 棒 性 能 好 等 特点 ， 
其 中 不 易 陷入 局 部 最 优 是 它 最 主要 的 优点 。 

本 文 主要 介绍 了 在 满足 系统 动态 响应 性 能 、 稳 定性 和 增益 
鲁 棒 性 的 基础 上 ， 利 用 频 域 法 获得 控制 器 的 公式 ， 并 作 图 求 出 
了 相应 控制 器 参数 ， 分 析 了 飞 蛾 火焰 优化 算法 的 特点 和 数学 模 
型 ， 详 细 阐 述 了 使 用 该 算法 优化 的 步 又， 进行 仿真 检验 ， 对 比 
和 分 析 了 两 种 方法 得 出 的 结果 。 最 后 ， 为 了 增加 实验 结果 的 可 
信 度 ,对 优化 后 的 控制 器 进行 方 波 和 正弦 波 速度 信号 跟踪 仿真 。 


1 “” 频 域 法 获得 控制 器 


相对 于 整数 阶 PI 控制 器 ， 分 数 阶 PT 控制 器 多 了 
参数 1， 能 在 更 大 范围 内 调节 控制 器 的 性 能 。 近 年 来 许多 学 者 
利用 分 数 阶 PI 控制 器 对 一 阶 和 二 阶 被 控 系统 进行 研究 ， 但 是 
现实 生活 中 的 被 控 系 统 实 际 上 就 是 一 个 分 数 阶 模型 ， 所 以 气 弃 
传统 的 整数 阶 被 控 系 统 ， 采 用 分 数 阶 被 控 系 统 更 为 合理 。 其 模 


个 可 池 


5 昌 


天 
IT (2) 
被 控 系 统 的 幅 值 和 相位 分 别 如 式 (3) (4) 所 示 。 
天 
P 7 le 
|P(j@®)| TY A + Bl) (3) 
| a | (4) 


其 中 : 


A(@) = o5 co CD] cos 和] (5) 


B(@) = ws sin CD sin [Ie] (6) 
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P(wO) =1+ 开 oO co 2 (9) 


Co) = Kw sin 4 (10) 


分 数 阶 PP 控制 器 的 设计 要 考虑 动态 响应 性 能 、 稳 定性 和 
增益 鲁 棒 性 的 要 求 ”"。 这 里 一 般 是 基于 以 下 三 个 指标 来 实现 
的 : 


a 增 益 穿越 频率 指标 。 根 据 系统 的 相 频 特性 和 幅 频 特 性 ， 
在 增益 穿越 频率 点 w 处 ， 开 环 传递 函数 的 幅 值 满足 


|cGw)PGw) 


=1 (1 


b) 相 位 裕 度 指标 。 为 了 保证 系统 的 稳定 性 ， 在 增益 截止 频 
率 点 下 的 相位 裕 度 9,, ， 其 满足 

Arg[CQJ@®.)PY®)]=-x+9, (12) 
c) 增 益 扰动 鲁 棒 性 指标 。 在 增益 截止 频率 @. 处 ， 系 统 的 相 
频 特性 曲线 的 导数 是 0。 此 时 意味 着 系统 对 外 界 的 干扰 具有 较 

强 的 鲁 棒 性 ， 具 体 的 表达 式 如 (13) 所 示 。 
dlArelCI%) PIO - 
40 w=0, 
当 给 定 了 系统 的 开 环 增益 穿越 频率 以 、 相 位 角 裕 度 9 、 控 
制 器 阶 次 4 ， 由 增益 和 鲁 棒 性 指标 可 得 


[ =0 (13) 


dL478(Cwo ,dlArePUo)] _, es 
dw Qn d@ =0, 
d[Arg(P(i 
引入 变量 M4 = 
0 = 
把 M 和 式 (8) 带 入 式 (14) 中 可 得 
一 4 天， mn[ 32]w 

+M=0 (15) 


wo K? +20 °K, cos [22] +1 


引入 变量 


N=2Mo@” co 32 一 asin( 4]jor” 


把 N 带 入 式 (15) 中 可 解 得 


(16) 
把 @.、44、; 带 入 式 (11) 可 以 解 得 


2 2 
TVA(@e) + B(@.) a 


” KVP(@.) +Q(@.) 


本 文采 用 的 被 控 对 象 是 PMSM， 具 体 的 电机 参数 如 表 1 所 


对 于 分 数 阶 PT 控制 器 式 子 式 (1) 的 幅 值 和 相位 如 下 所 示 : 
[CUO K,NP(O) +Q(0) (7) 


es _C(w) 
47s[C(Jw)] = sen| 2 | 


(8) 


其 中 : 


示 。 利 用 上 述 的 电机 参数 ， 采 用 数值 拟 合 的 方法 对 其 电磁 环 和 
机 械 环 进行 分 数 阶 系 统 时 域 参数 辨识 ， 可 以 求 得 被 控 对 象 的 分 
数 阶 模型 为 


6196.69 
S15 十 325.298108 +3974.66 


式 (18) 可 以 得 到 被 控 对 象 的 截止 频率 @.=12 rad/s， 相 位 


G(s)= (18) 
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角 裕 度 9 =90" ， 时 间 常 数 六 0.8 s。 
分 别 画 出 式 (14)(15) 的 图 形 。 关于 4 的 函数 图 像 ， 如 图 1 
所 示 。 
表 1 永 磁 同 步 电 机 的 参数 
电机 参数 参数 值 
额定 转速 /r/min) 2500 
定子 电阻 RM 3.987 
定子 dg 轴 电 感 LH 0.0088 
转动 惯量 (ke*m?) 0.0075 
励磁 磁 通 Wb 0.175 
极 对 数 p 1 


频 域 法 获得 控制 器 参数 图 


分 数 阶 积 分 参数 


0.6 0.7 0.8 0.9 1 1.1 2 1.3 1.4 1.5 
分 数 阶 阶 次 lamda 


图 1 


频 域 法 获得 控制 器 的 


图 中 红线 是 式 (14) 的 波形 图 ， 蓝 线 是 式 (15) 的 波形 图 


A 
， 巴 


们 的 交点 为 所 求 的 参数 ， 则 (1 ， 
(17) 可 以 得 到 必 =3.672， 进 而 得 到 分 数 阶 PT 控制 器 为 


C(s) = 3.6720 + 2 


五) 的 值 为 (0.857,12.22)， 带 入 式 


(19) 


利用 系统 被 控 的 分 数 阶 模型 式 (18) 和 频 域 法 得 到 的 PT 控 


一 上 


所 器 式 (19), 可 以 得 到 系统 的 闭环 传递 函数 , 对 


其 闭环 传递 函 


数 


进行 阶 跃 仿真 , 则 可 以 得 到 系统 的 阶 跃 响应 波形 ,如 图 


Step response 
T 


i 1 1 
06 08 1 12 14 16 18 2 
Time [sl 


图 2 频 域 法 得 到 的 系统 阶 跃 响应 波形 


1 1 
0 02 04 


从 图 2 


的 阶 跃 响应 波形 


5 


2 所 示 。 


可 以 得 到 ， 上 升 的 时 间 刀 =0.12 s， 
骨节 的 时 间 大 约 为 大 =0.6s， 超 调 量 为 4.76%。 虽 然 此 时 得 到 的 


分 数 阶 控 


的 增益 
应 跟随 
大 
止 频率 
索 最 


多 的 增益 穿越 频率 %. 和 分 数 阶 次 4 。 
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一 一 


制 器 可 以 确保 特定 的 增益 穿越 频率 、 相 位 裕 度 和 不 错 
扰动 鲁 棒 性 ， 但 是 仍 不 能 使 控制 系统 达到 最 佳 的 阶 跃 响 
性 能 

此 分 数 阶 控制 器 设计 就 转换 成 了 如 何 得 到 最 优 的 增益 截 
@. 和 分 数 阶 阶 次 4 的 问题 。 飞 蛾 火焰 优化 算法 被 用 来 搜 
将 每 一 次 得 到 的 


(w,W) 代入 公式 求 得 分 数 阶 速度 控制 器 ， 并 在 MATLAB 中 进 


行 阶 跃 
计 控 制 | 
值 最 小 


响应 仿真 。 同 时 时 间 乘 误差 平方 积分 (ITAE) 被 计算 来 估 
器 的 性 能 。 通过 飞 蛾 火焰 优化 算法 的 迭代 计算 , 当 Jras 
时 ， 得 到 的 分 数 阶 控制 器 是 最 佳 的 。 


2 “ 飞 蛾 火焰 优化 算法 (MFO) 


2.1 

不 
优化 过 
焰 优 
机 制 的 


飞 蛾 火焰 优化 算法 的 特点 


同 结构 的 控制 器 ， 参 数 优化 的 方法 也 有 所 不 同 。 为 了 使 
程 不 易 陷入 局 部 最 优 ， 并 提高 优化 效率 ， 本 文 利用 飞 蛾 
化 算法 来 优化 分 数 阶 控 制 器 ， 该 算法 是 受 飞 蛾 定位 导航 


启发 ， 通 过 飞 蛾 火焰 螺旋 飞行 数学 模型 更 新 飞 蛾 火焰 位 


优化 过 


a) 螺 旋 的 初始 点 从 飞 蛾 初始 化 空 


空间 位 


MFO 算法 选择 螺旋 函数 作为 飞 峨 空间 


最 终 收敛 于 火炮 位置 的 新 型 群 智能 仿生 算法 。 
立 置 的 更 


新 算 子 , 在 


程 中 遵循 以 下 准则 : 


则 位 置 开始 ， 火 焰 的 最 佳 


是 螺旋 的 最 后 一 点 。 


b) 螺 旋 的 波动 范围 不 超出 搜索 范围 
co) 在 火焰 周围 ， 


而 收敛 


& 蛾 能 够 随 tt 为 [-1,1] 间 的 随机 数 ) 的 变化 
到 任何 点 ,，t 越 小 ， 离 火焰 的 距离 越 近 


dd) 当 飞 蛾 靠近 火焰 时 ， 火 焰 两 侧 的 位 置 更 新 频率 增加 。 


2.2 飞 蛾 火焰 优化 算法 的 数学 模型 
MFO 算法 数学 模型 描述 如 下 所 示 : 


mi mj; mja 
11l27 Mm,, mya 
M= : 1 1 (20) 
m1 mm, » Ma 
OM=[OM, OM, :… OH] (21) 
在 矩阵 式 (20) 中 ,n 为 飞 蛾 的 种 群 规模 ;，4d 为 优化 问题 维度 
的 飞 蛾 所 处 的 空间 位 置 。 和 矩阵 式 (21) 是 用 来 存储 飞 蛾 个 体 的 适 
应 度 值 。 
MFO 算法 中 另 一 个 关键 部 分 是 火焰 ， 表 示 为 
FF, Fi, Fi | 
EF, EF,, Fy 
F=| *” 2 (22) 
Fj FF, Fa 过 
并 利用 式 (23) 和 矩阵 存储 火焰 适应 度 值 ， 即 
OF =[OF, OF, :. OF,| (23) 


其 中 : 


n 为 飞 蛾 种 群 规 模 ;，4 为 优化 问题 维度 。 
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飞 蛾 和 火炮 都 是 候选 解 ， 它 们 的 主要 区 别 是 在 进化 过 程 中 
位 置 的 更 新 方式 不 同 。 飞 蛾 是 在 搜索 空间 附近 移动 的 搜索 动因 ， 
火焰 是 当前 迭代 所 获得 的 最 佳 位 置 。 通 过 火焰 位 置 不 断 地 更 新 
直至 获得 最 佳 的 位 置 。 因 此 ，MFO 算法 可 以 定义 为 

MFO=(1,P,K) (24) 


i t 


其 中 :了 为 随机 生成 的 飞 蛾 种 群 规模 和 相应 的 适应 值 函 数值 ; 
P 为 飞 蛾 在 搜索 空间 移动 的 函数 , 其 为 嘱 和 矩阵 和 位 置 更 新 返回 值 ; 
KK 为 终止 条 件 是 否 满足 的 判别 函数 。 
用 数学 模型 描述 飞 蛾 火焰 位 置 更 新 如 下 : 
M,=S(M,,F,) (25) 


其 中 : M; 表 示 第 i 只 飞 峨 位 置 ， 了 表示 第 j 个 火焰 !s 为 螺旋 


定义 用 于 模拟 飞 蛾 螺旋 飞行 路 径 的 MFO 算法 表达 式 如 下 : 


S(M,,F,)= Die” cos(2xt) +F, (26) 


式 (26) 是 对 飞 蛾 螺旋 曲线 盘旋 运动 的 路 径 模 拟 ， 
迹 如 图 3 所 示 。 


其 几何 轨 


0.5 -0.5 = 站 0 0.5 耳 


图 3 ” 飞 蛾 绕 对 数 螺旋 曲线 几何 轨迹 图 


其 中 : S 表示 飞 蛾 Mi 绕 火 焰 五 运动 更 新 后 的 位 置 ， 
DF 一 Mi| 是 第 i 只 飞 蛾 MM; 到 第 i 个 火焰 五 的 距离 ;5b 为 螺 
旋 形 状 常数 ;z 为 [-1.1] 之 间 的 随机 数 ， 它 定义 了 飞 蛾 距离 火焰 
下 一 位 置 的 远近 程度 。 

为 保证 MFO 算法 获得 较 快 的 收敛 速度 ， 提 出 自 适 应 火焰 
IE 


N-1 
赤 ) (27) 


flameno =round(N 一 1 


其 中 : /为 当前 迭代 次 数 ; N 为 最 大 火焰 数量 ; 
数 。 
2.3 ” 飞 蛾 火焰 优化 算法 的 设计 和 流程 

选取 ITAE 作为 系统 动态 性 能 的 指标 , 搜索 使 ITAE 值 达 到 
最 小 值 的 控制 器 。 ITAE 指标 具有 很 好 的 选择 性 ， 系统 参数 的 微 
小 变化 能 够 显著 改变 ITAE 值 。ITAE 指标 如 式 (28) 所 示 。 


T 为 最 大 达 代 次 


J rg = | | el(1) kt (28) 
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其 中 : e( 的 为 系统 期 望 输出 与 实际 输出 之 间 的 偏差 。 


设计 的 算法 流程 如 图 4 所 示 。 


用 频 域 法 求 出 寻 优 范围 
的 中 心 值 


确定 待 优化 参数 的 维 
度 、 搜 索 范 围 和 算法 迁 
代 终 止 条 件 


Y 
确定 MF0 算 法 的 适应 度 
函数 


v 


计算 飞 蛾 个 体 适应 度 函 
数值 ， 并 找到 当前 最 好 
飞 蛾 个 体位 置 并 将 其 保 
存 为 火焰 适应 度 值 矩阵 


判断 算法 迭代 终止 条 件 


是 否 满足 
| ] 
是 
迭代 开始 ， 更 新 火焰 数 


量 ， 计 算 火 烙 与 飞 峨 间 
的 距离 ， 并 更 新 飞 蛾 火 
焰 位 置 


一 


i 
计算 飞 蛾 个 体 适 应 度 函 最 终 
类 的 空间 位 置 及 所 对 应 的 
0 St 
YY | 
投 到 当前 最 好 上 如 个 体 算法 结束 
4 算法 流程 
飞 蛾 火焰 优化 算法 的 具体 参数 设置 如 表 2 所 示 。 
表 2 飞 峨 火焰 优化 算法 参数 
参数 名 参数 值 
飞 峨 种 群 规 模 100 
最 大 火焰 数量 N 100 
最 大 迭代 次 数 500 
螺 线 塑 形 常 数 1.5 
终止 条 件 己 知 的 最 优 解 或 者 评价 的 次 数 超过 500 次 


3 ”仿真 检验 


3.1 求解 控制 器 
永 磁 同步 电机 分 数 阶 速度 控 
模型 。 


制 系 统 可 以 为 图 5 所 示 的 简化 


图 5 永 磁 同步 电机 控制 系统 分 数 阶 简化 模型 
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图 中 ，P(s) 为 永 磁 同步 电机 速度 控制 系统 的 被 控 对 象 ， 为 
了 使 仿真 更 加 的 合理 ， 被 控 对 象 采用 分 数 阶 模型 ，C(s) 为 分 数 
阶 的 速度 PE 控制 器 

采用 前 面 提 到 的 频 域 法 和 飞 蛾 火焰 优化 算法 对 分 数 阶 PE 
叫嚣 参数 进行 整定 ， 其 中 飞 蛾 火焰 算法 的 迭代 过 程 如 图 6 所 


i 


Convergence curve 
- - 


= T T 
10 


Best flame (score) obtained so far 


100 200 300 A400 S00 
lteration 


收 火 焰 优 化 算法 的 迁 代 过 程 


图 6 飞 


为 了 模拟 飞 蛾 火焰 优化 算法 的 迭代 过 程 和 和 
焰 优 化 算法 进行 测试 函数 的 仿真 ， 如 图 7 所 示 。 


吉 果 ， 对 飞 蛾 火 


Test function 


到 


7 怀 蛾 火焰 优化 算法 的 测试 函数 仿真 结果 
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图 8 火焰 优化 算法 优化 的 分 数 阶 控制 系统 模型 


可 以 得 到 控制 系统 的 bode 图 ， 如 图 9 


式 (18) 和 (29) 念 


所 示 。 


Bode Diagram 


a 
S 
吕 

革 


己 
3 


Frequency (rad/s): 34.4 


System: H1 am | 
Magnitude (dB): -0.106 | 
| 


Magnitude (dB) 
o 


豆 
3 


-200 5 ! 


四 | 
System H1 
Frequency (rad/s): 34.4 
Phase (deg): -105 中 


Phase (deg) 


10° 10" 10° 
Freauency (rad/s) 


图 9 采用 C(s) 的 控制 系统 的 波 特 图 


将 图 9 的 幅 频 特性 曲线 和 相 频 特性 曲线 的 参数 进行 整理 ， 
如 表 3 所 示 。 
表 3 PT 控制 系统 的 bode 图 参数 
角速度 rad/s 相位 deg 幅 值 dB 线 角 /* 
34.4 -105 -0.106 0 
从 图 9 和 表 3 可 以 看 出 ， 系 统 的 相 频 特性 曲线 在 增益 频率 


点 处 是 水 平 的 。 所 以 ， 当 系统 增益 出 现 微小 的 扰动 时 ， 系 统 的 
相位 裕 度 能 基本 保持 不 变 。 也 就 是 说 ， 控 制 系统 的 相位 裕 度 满 
足 边 界 条 件 的 要 求 ， 控 制 系统 具有 一 定 的 鲁 棒 性 。 
3.2 ”优化 后 的 分 数 阶 控制 器 与 频 域 法 控制 器 比较 

优化 后 的 分 数 阶 控制 器 与 频 域 法 控制 器 比较 如 图 10 所 示 。 
把 两 个 系统 的 阶 跃 响应 的 性 能 指标 以 表格 的 形式 进行 对 比分 析 ， 


通过 图 6 可 以 看 出 ， 随 着 迭代 次 数 的 增加 ， 目 标 函 数 的 适 
应 度 值 不 断 地 减 小， 说 明 优化 的 过 程 是 比较 合理 的 。 在 迭代 次 
数 到 达 180 次 左右 时 , 目标 函数 基本 达到 了 稳定 状态 .迭代 180 
次 时 已 经 确定 了 最 优 值 ， 同 时 相对 于 给 定 的 迭代 次 数 500 次 ， 
其 途 代 结 果 是 理想 的 。 

通过 飞 蛾 火焰 算法 优化 后 得 到 @. =34.4， 
=15.58， 带 入 式 (17) 可 以 得 到 上, =2.003; 则 分 数 阶 PT 控 
式 (29) 所 示 。 


4=1.023。 则 kk 
制 器 如 


31.2067 


C(s) =2.003+ rm (29) 


根据 式 (18) 和 (29) 可 以 建立 分 数 阶 PE 控 制 
算法 的 优化 仿真 模型 ， 如 图 8 所 示 。 


器 飞 蛾 火焰 优化 


如 表 4 所 示 。 


Step response 
T T T 


;一 一 顶 花 后 的 SI 


Amplitude 


1 1 h h h 
05 06 07 08 09 1 


Time [s] 


1 i h 
0 01 02 03 04 


图 10 频 域 法 与 飞 蛾 火焰 算法 优化 比较 
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表 4 两 种 控制 器 参数 对 比 不 。 
控制 器 。 ”调节 时 间 L/s ”上 升 时 间 /s 。” 超 调 量 o% 表 5 控制 器 的 阶 跃 响应 性 能 指标 
G(1) 0.3 0.11 4.60 控制 器 峰值 / 转 速度 稳定 时 刻 /s ITAE 值 
G2) 0.5 0.15 4.32 C(s) 1060 0.25 1.038 
Ci(s) 1090 0.16 4.567 
从 图 10 和 表 4 可 以 看 出 ， 优 化 得 到 的 分 数 阶 PT 控制 器 使 CC DG 1 Gs 
系统 的 调节 时 间 和 上 升 时 间 明 显 降低 ， 从 而 提高 了 系统 的 动态 CC 0 Ge O83 
响应 性 能 。 从 超 调 量 这 一 方面 来 看 ， 优 化 的 控制 器 与 频 域 法 相 CGO) 1070 1 
差 不 大 , 以 后 可 以 通过 改进 飞 峨 火焰 算法 来 提高 这 方面 的 能 力 。 
3.3 优化 的 分 数 阶 控制 器 之 间 的 比较 从 图 11 和 表 5 可 以 看 出 ， Cs) 的 响应 速度 比 Cx(s) 慢 ， 
于 控制 系统 的 相位 裕 度 已 经 满足 边界 条 件 的 要 求 了 ， 所 ”但 是 超 调 量 、 波 动 和 ITAE 相对 较 小 ， 从 而 可 以 得 出 结论 C1(5) 
以 为 了 验证 获得 的 控制 器 为 最 优 的 .从 图 9 中 的 最 佳 增益 频率 ” ”控制 器 的 阶 跃 响应 性 能 要 优 于 C(s) 。 同 理 ， 也 可 以 看 出 C3(s) 
和 相位 裕 度 的 附近 选择 四 组 分 数 阶 PT 控制 器 ， 并 将 这 五 组 优 ”的 响应 速度 比 C4(s) 快 ， 并 且 超 调 量 、 波 动 和 ITAE 更 小 ， 因 此 
化 的 控制 器 进行 电机 速度 跟踪 仿真 测试 比较 。 可 以 证 明 以 上 的 验证 方案 是 合理 的 。 
分 数 阶 控制 器 C(s) 采用 前 面 求 得 的 增益 截止 频率 从 图 11 和 表 5 可 以 看 出 ， 尽 管 CCs) 和 Cs(s) 控制 系统 的 
@.=34.4rad 1s 和 相位 裕 度 ov =75 ,在 其 附近 选取 的 四 组 。 输出 能 够 更 快 地 达到 稳定 状态 ， 但 是 超 调 量 和 ITAE 都 比 C(s) 
(@.,9,) 分 别 为 大 , 而 Cs(s) 和 C4(s) 控制 系统 的 输出 响应 速度 要 比 C(s) 慢 , 且 
C1l: @&.=40rad/s 70 超 调 量 和 ITAE 都 比 C(s) 大 , 因此 使 用 所 得 到 的 C(s) 控制 系统 
C2: @.=45rad/s =65 拥有 更 为 优异 的 动态 响应 性 能 。 
C3: 0@®.=30rad/s 80 3.4 鲁 棒 性 验证 
C4: 中 .=25rad/s 9, =85 为 了 验证 优化 后 的 分 数 阶 系统 具有 较 强 的 鲁 棒 性 ， 将 分 数 


很 明显 ,C1 比 C2 更 接近 


基于 这 四 组 值 ， 可 以 求 出 对 应 的 四 组 分 数 阶 控制 器 


33.8890 


so 932 


C(s)=2.164+— oO— 


37.8529 


eT 


C,(s)=2.234+ 


30.7500 


1.109 
RY 


C,(s) =1.885+ 


28.5339 


st 195 


Ci(s)=1.676+ 


、C2(S) 、C3(5) 、G4(5) 来 控 


记 电 机 的 速度 为 1000 wmin, 在 MATLAB 


上 分 别 用 


(Q@.,9,) ,C3 比 C4 更 接近 (@.,9,)。 


制 分 数 阶 永 磁 同步 


Amplitude 
只 


线 。 仿 真 的 结果 如 图 


它们 的 速度 阶 跃 响应 


Step response 


11 所 示 。 


1200 F F F F F TF F TF F 


1000 /77 


了 下 了 让 了 了 
0 01 02 03 04 05 06 07 08 09 
Time [s] 


图 11 


图 11 的 不 同 控制 器 的 阶 跃 响应 性 能 指标 列强 


不 同 分 数 阶 控制 系统 的 速度 跟踪 响应 仿真 


,如 表 5 所 


C(s) 、 


阶 PP 控制 器 的 比例 系数 分 别 设置 为 最 
80%， 即 2.003 0、2.403 6、 
包机 模型 进行 速度 阶 跃 响应 仿真 测试 。 


佳 优化 值 的 100%、120%、 


1.602 4， 并 将 三 组 控制 器 对 分 数 阶 


三 组 控制 器 为 
C=2.0030 +3 .2007 
ed 
ns 
C=1.6024 + 30 


仿真 结果 如 图 12 所 示 。 


Step response 
Tt T 


Ss 


1200r F Fr F r 十 
1000 上 一 | Kp=2.0030 | 
Kp=2.4036 
kp=1.6024 
800 上 - [一 kp=16024 | 
0 
百 
己 
考 600- 
E 
< 
400- + 
200|- 
0 005 01 015 02 025 03 035 04 045 05 
Time [s] 
图 12 优化 的 分 数 阶 控制 器 鲁 棒 性 验证 


制 器 鲁 棒 性 性 能 指标 如 表 6 所 示 。 
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表 6 控制 器 鲁 棒 性 性 能 指标 如 图 14 所 示 ， 将 幅 值 为 42r /min ， 频 率 为 20 Hz 的 正弦 
Ks 什 超 调 量 /% 速度 稳定 时 刻 /s 童 号 作为 控制 系统 的 参考 速度 信号 ， 将 其 作用 在 优化 后 的 分 数 
2 0030 46 0.255 阶 控制 系统 中 。 其 中 速度 是 时 刻 在 变化 的 ， 没 有 稳 态 的 时 刻 ， 
2 4036 4.63 0.255 因而 波形 呈现 一 种 滞 环 的 形式 ， 这 是 控制 器 动态 调节 的 结果 。 
1.6024 4.58 0.255 超 调 量 约 为 1%， 这 是 因为 正弦 波 两 个 相 邻 时 刻 的 速度 相差 不 
大 ， 所 以 调节 过 程 相 比 阶 跃 信号 比较 快 。 
从 图 12 和 表 6 可 以 看 出 , 三 条 曲线 的 超 调 量 很 接近 , 同时 所 以 从 超 调 量 和 整体 的 波形 来 看 ， 优 化 的 控制 系统 能 够 较 
速度 的 稳定 时 刻 也 基本 相同 ， 因 此 可 以 得 出 结论 ， 增 益 扰 动 大 好 地 跟踪 正弦 速度 信号 。 
约 在 最 佳 值 的 +20% 时 ,控制 系统 的 相位 裕 度 也 几乎 是 相同 的 ， 4 ”结束 语 
这 足够 满足 系统 对 增益 扰动 的 鲁 棒 性 要 求 。 
3.5 ”优化 后 的 控制 器 方 波 和 正弦 波 速度 信和 号 跟踪 仿真 本 文 提出 了 一 种 飞 峨 火焰 算法 优化 分 数 阶 控制 器 的 方法 。 
前 面 已 经 验证 了 控制 系统 具备 良好 的 动态 响应 性 能 ， 为 了 对 于 被 控 对 象 ， 据 弃 传 统 的 一 阶 和 二 阶 模型 ， 采 用 更 为 合理 的 
进 人 分 数 阶 模型 作为 被 控 对 象 。 首 先 使 用 频 域 法 得 到 分 数 阶 PP 控 
速度 设置 为 方 波 信 号 和 正弦 波 信号 。 事 实 上 ， 方 波 信号 和 正弦 制 器 的 参数 ， 然 后 用 算法 对 该 分 数 阶 PP 控制 器 参数 进行 优化 。 
波 信号 可 以 看 成 无 穷 个 阶 跃 信 号 变化 的 过 程 。 具 体 的 波形 如 图 为 了 验证 算法 得 到 的 控制 器 的 优越 性 ， 将 其 与 频 域 法 得 到 的 控 
13、14 所 示 。 制 器 进行 仿真 对 比 ; 为 了 验证 最 优 控 制 器 的 动态 响应 性 能 ， 将 
方 流速 度 信号 中路 其 与 周围 的 值 形成 的 控制 器 进行 对 比 验证 ;为 了 验证 其 增益 鲁 
L i 棒 性 ， 通 过 对 比 k; 在 +10% 的 范围 变化 时 的 响应 曲线 ， 最 后 对 
~ 方 波 和 正弦 波 速 度 信号 进行 仿真 。 仿 真 结果 表明 ， 优 化 后 的 分 
i 数 阶 系统 具有 很 好 的 动态 响应 性 能 、 较 强 的 和 鲁 棒 性 和 不 错 的 稳 
归 。 定性 ; 但 是 优化 的 控制 器 与 频 域 法 相 比 ， 超 调 量 并 没有 很 大 的 
站 -5 改善 ， 这 也 有 待 以 后 对 飞 蛾 火焰 算法 在 捕 焰 、 弃 焰 以 及 适应 度 
号- 函数 等 方面 的 改进 。 
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